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Abstract 

We present a preliminary measurement of the CP-violating parameters in fully reconstructed 
B'^^D^*^^tt^ and B'^^D^ decays in approximately 232 million T(4S') BB decays collected 
with the BABML detector at the PEP-II asymmetric-energy B factory at SLAC. From a maximum 
likelihood fit to the time-dependent decay distributions we obtain for the parameters related to the 
CP violation quantity sin(2/3 -|- 7): 

^Dtt ^ -0.013 ±0.022 (stat.) ±0.007 (syst.) , cgp = -0.043 ± 0.042 (stat.) ± 0.011 (syst.) 
^D-TT ^ -0.043 ± 0.023 (stat.) ± 0.010 (syst.) , cg*'^ = 0.047 ± 0.042 (stat.) ± 0.015 (syst.) 
a'^p = -0.024 ± 0.031 (stat.) ± 0.010 (syst.) , = -0.098 ± 0.055 (stat.) ± 0.019 (syst.) 
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1 INTRODUCTION 



While the measurement of sin 2/3 is now a precision measurement 012]) the constraints on the other 
two angles of the Unitarity Triangle a and 7, are still limited by statistics and by theoretical 
uncertainties. We report on the measurement of CP-violating asymmetries in B^^D^*^^Tr^ and 
B^^D^p^ decays in T{4S) —>■ BE decays, which are related to sin(2/3 + 7) j^Ej. This analysis 
updates results already presented in Ref. [HI using a larger data sample (232 million instead of 110 
million T{AS) — > BE decays), and complements the recently published BABAR measurement of CP 
violation in i? — > D*^tt^ decays with a partial reconstruction technique [Jj. 

The time evolution of B^^D^*^^h^ decays, where is a charged meson made of u and d 
quarks, is sensitive to 7 because the CKM-favored decay amplitude B^^D^*^^h~ , which is pro- 
portional to the CKM matrix elements V^^^V*^, and the doubly-CKM-suppressed decay amplitude 
B^^D^*^~^h~ , which is proportional to VcdV*^^, interfere due to B'^B^ mixing. The weak phase 
difference between the two decay amplitudes is 7. When combined with B^B^ mixing, the total 
weak phase difference between the interfering amplitudes is 2/3 -|- 7. 

The decay rate distribution for B^^D^*^^h^ decays, neglecting the decay width difference, is 
given by 

/±(?7,At) = x[lTS^sm{AmdAt)^riCcos{AmdAt)], (1) 

where r is the B^ lifetime, Am^ is the B^B^ mixing frequency, and At = tree — ^tag is the time of 
the B'^ — > D^*^^h^ decay (i?rec) relative to the decay of the other B (i?tag)- In this equation the 
upper (lower) sign refers to the flavor of i?tag as B^ (-B"), while rj = +1 (—1) and C = + (~) for the 
final state D^*^^h'^ {D^*^^h^). In the Standard Model, the S and C parameters can be expressed 
as 
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2Im(A±) , ^ 1-^ 

S± = - , , |\ 10 , and C=-—^, (2) 



1 + |A±|2' 1 + r 



where r = |A+| = 1/|A_| and 



pAiBO^DTh^)-"^ ■ 

Here ^ is a function of the elements of the mixing Hamiltonian 8^, and 6 is the relative strong 
phase between the two contributing amplitudes. In these equations, the parameters r and 5 depend 
on the choice of the final state and will be indicated as r^'^, J^'", in the B^^D^TT^ case, and 
r^P, 6^^, in the B^^D^p^ case. For the B^^D*^Tr^ mode the expression is similar with r^*'^ 
and 6^*'^ (according to Ref. [S| the strong phase is 6^*"^ -|- vr in this mode, but this does not 
affect this measurement). Interpreting the S and C parameters in terms of the angles of the 
Unitarity Triangle requires knowledge of the r parameters. The r parameters are expected to be 
small (~ 0.02) and therefore cannot be extracted from the measurement of C with the current 
statistics. They can be estimated, assuming SU{3) symmetry and neglecting contributions from 
annihilation diagrams, from the ratios of branching fractions 13{B'^—)-Di*^^Tr~)/B{B^—i-D^*'^~Tr~^) 
and B{B^^D+p-)/B{B^^D'p+) dEUT]. Note that the B°^D*+tt~ and B^^D+p~ decays [S] 
have not yet been observed. 
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2 THE BABAR DETECTOR AND DATASET 

This measurement is based on 232 million T(4S') BB decays collected with the BABAR detector 
at the PEP-II storage ring. The BABAR detector is described elsewhere |1U| . We use Monte Carlo 
simulation of the BABAR detector based on GEANT4 JJ| to validate the analysis procedure and 
to estimate some of the backgrounds. 

3 ANALYSIS METHOD 

The analysis strategy is similar to our previous publications on this topic [HI El and to other 
time-dependent measurements performed at BABMl pp. 

We reconstruct the B^^c in the D^vr^, D*^tt^ and D^p^ final states. Candidate B D*^tt^ 

decays are reconstructed with the D*^ decaying to i)^7r^, where the subsequently decays to 
one of the four modes K^-k^, K^-k^tt^ , K^'k^tt^'k^ , or K^^-K^-K^. Candidate B D^tt^ and 
B — > D^p^ decays are reconstructed with the decaying into K^tt^tt^ and K^ir^, and the 

decaying into vr^vr'^. The D^*^^ candidates are then combined with either a single track or a 
track and a vr candidate with invariant mass in the window, 620 < m(7r±7r°) < 920 MeV/c^. 
Exploiting the spin properties of the decay of a pseudo-scalar meson into a pseudo-scalar and a 
vector meson, we require the cosine of the helicity angle defined as the angle between the 
charged pion and the D momentum in the p^ rest frame, to satisfy |cos0hei| > 0.4. No requirement 
is applied on any helicity angle of the B — > Z'*^7r^ decay channel. The event selection and candidate 
reconstruction are described in more details in Refs. pi I12|. 

To identify the flavor of the Btag, we use multivariate algorithms that identify signatures in the 
B decay products that determine ("tag") the flavor to be either a B^ or a B^ . Primary leptons from 
semi-leptonic B decays are selected from identified electrons and muons and from isolated energetic 
tracks. The charges of identified kaons and soft pions from D*^ decays are also used to extract flavor 
information. These algorithms are combined, taking into account the correlations among different 
sources of flavor information, and provide an estimate of the mistag probability. The tagging 
procedure has been improved with respect to the procedure used in our previous analysis [H] with 
the addition of information from low momentum electrons, A ^ pvr decays, and correlations among 
identified kaon candidates. Each event with mistag probability less than 45% is assigned to one 
of six hierarchical, mutually exclusive tagging categories. The lepton tagging category contains 
events with an identified lepton, while other events are divided into categories based on the mistag 
probability. The effective efficiency of the tagging algorithm, defined as Q = ^itii). — 2u)j)^, where 
ej and Wi are the efficiency and the mistag probability for the tagging category i, respectively, 
improves by 5% (relative) over the algorithm used in Ref. [5]. 

The time difference At is calculated from the measured separation along the beam collision axis 
Az between the Brec and Btag decay vertices. We determine the B^^c vertex from its charged tracks. 
The Btag decay vertex is obtained by fitting tracks that do not belong to Brec, imposing constraints 
on the Brec momentum and the beam-spot location. The At resolution is approximately 1.1 ps. 

Signal and background are discriminated by two kinematic variables: the beam-energy sub- 
stituted mass, mES = \J (\/^/2)^ — P*b^ ^ and the difference between the B candidate's measured 
energy and the beam energy, /S.E = E'^ — (y^/2). E'^ the energy (momentum) of the B 

candidate in the e^e~ center-of-mass frame, and ^/s is the total center-of-mass energy. The signal 
region is defined to be |A£'| < 3a, where the resolution a is mode-dependent and is approximately 
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20MeV as determined from data. Figure Q shows the m^s distribution for candidates in the signal 
region. 

The background remaining after the event selection can be separated into two categories, one 
of which is due to random combinations of particles in the event {combinatorial background) and 
the other is due to B decays into final states similar to the signal final states, which have an rriES 
distribution similar to the signal {peaking background). 

The niES distribution of signal events is described by a Gaussian having a width of about 
2.5MeV/c^, whereas the peaking background is described by a bifurcated Gaussian with widths of 
about SMeV/c^ on the left side and 2.5MeV/c^ on the right side. To separate the combinatorial 
background, the mES distribution is fitted with the sum of a threshold function and the signal 
and peaking background components. 

We estimate the contributions of background peaking in the mss signal region using simulated 
events. The most dominant peaking background modes are B decays into open-charm final states 
similar to that of the signal, e.g. B^ D*^K~^, B^ Z)*°7r+//9+ or B^ D*~ p'^ . The Gaussian 
yields and the amounts of peaking background are summarized in Table ^ The latter are identified 
by their source as either coming from neutral or charged B meson decays. In the case of the 

Table 1: Signal yields, purities P, and fractions of peaking backgrounds fpeak of the selected samples 
for events with tagging information. The purity is computed in the signal region, for nriEs in a 



three a region around the B^ mass nominal va 


ue. 




Decay 


yields 


P{%) 


fpeak{%) 




mode 






B^ 






B^D^TT^ 


15635±135 


85.5±0.3 


3.3±0.1 


1.3±0.1 




B^D*^TT^ 


14554±126 


93.0±0.2 


2.8±0.1 


0.8±0.1 




B^D^p^ 


8736±105 


81.7±0.4 


1.3±0.2 


1.5±0.2 



B^ D p^ decays, an additional source of background must be considered, which has the same 
final state, B^ — > D^tt^tt^, where the tt^tt^ system is not produced through the p resonance. This 
background can introduce a dependence of the A parameter of Eq. |31on the vrvr'^ invariant mass: its 
contribution has been studied in Ref. [5]. This part of the analysis is not updated in this paper; 
we use our old result in the evaluation of the systematic error (section ^ . 

An unbinned maximum likelihood fit is performed to the time distribution of events in this 
sample. Taking into account possible CP violation effects on the tag side the likelihood for 
signal events of Eq. 1 for each tagging category i and for each decay mode = D^tt^ , D*^Tr^ , D^p^ 
becomes: 

/f '^(r/. At) = — X [1 ^ (a'^ ^ i]bi - ri4) sin(AmrfAt) ^ r? cos(AmdAt)] , (4) 

4r 

where in the Standard Model 

a^' = 2r'' sin(2/3 + 7) cos J'' , 
hi = 2r • sin(2/3 + 7) cos (5- , 

4 = 2cos(2/3 + 7)(r^sin5^-r^sin(5-). (5) 

Here r[ {5[) is, for each tagging category, the effective amplitude (phase) used to parameterize the 
tag side interference. Terms of order r^^ and r'^ have been neglected. Results are quoted only 
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7000 p 




Figure 1: tties distributions in the signal region for, from top to bottom, the B D tt^, 
B iD*^7r^ and B D^p^ sample for events with tagging information. Each solid curve repre- 
sents the result of a fit with a Gaussian (that describes signal events) plus a bifurcated Gaussian 
(that describes peaking background events) plus a threshold function (that describes combinatorial 
background events). Each dashed curve represents the background events only. 
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for the six and cj^^ parameters, which are independent of the unknown parameters r'- and 5'^ 
(semi-leptonic B decays have no doubly CKM-suppressed contributions and therefore r^'^^ = 0). 
The other parameters are constrained by the fit, but, as they depend on r'- and 5-, they do not 
contribute to the interpretation of the result in terms of sin(2/3 + 7). 

The signal At distribution in Eq.^is convoluted with the resolution function parametrized with 
the sum of three Gaussians to take into account finite At resolution. The probabilities of incorrect 
tagging (wi) are accounted for by multiplying the o'^, cf parameters and the cos(AmrfAt) term 
by the dilutions Di = I — 2wi (average dilutions for B^ and B^). Possible differences in mistag 
probabilities for B^ and B^ are also taken into account. The resolution function and the tagging 
parameters are floated in the fit to the data and are consistent within errors with previous BABMi 
analyses p. 

The likelihood function has a contribution for each source of background. The combinatorial 
background is parametrized as the sum of a component with zero lifetime and one with an effective 
lifetime floated in the fit to the data. The fraction of each background component is determined 
from the events in the tties sidebands, having 5.20 < rriEs < 5.27GeV/c^. The background events 
are described using effective dilution parameters, obtained from the fit to the data, and have a At 
resolution function consisting of the sum of two Gaussians, also floated in the fit to the data. The 
charmed peaking background coming from B^ mesons is modeled by an exponential with the B^ 
lifetime, and its amount is fixed to the value predicted by simulation. The resolution function is the 
same as the signal resolution, while the dilution parameters are fixed to the values obtained from 
a B~^ control sample. The charmed and charmless peaking backgrounds from B^ mesons, whose 
amounts are also fixed to the value estimated from simulation, are modeled with a likelihood similar 
to the signal likelihood with no CP violation (all the a, 6, c parameters fixed to 0). Possible CP 
violation in this background is taken into account in the evaluation of the systematic uncertainties. 
The resolution and the dilution parameters are the same as for the signal. 

4 SYSTEMATIC STUDIES 

Table 121 shows the contributions to the systematic uncertainties of the a and Qep CP parameters. 



Table 2: Systematic uncertainties on the a and c parameters. 







B'^^D*^TT^ 




Source 














Vertexing (aAt) 


0.0037 


0.0064 


0.0080 


0.0110 


0.0047 


0.0110 


Fit (cjfit) 


0.0051 


0.0088 


0.0052 


0.0093 


0.0075 


0.0129 


Model (fJmod) 


0.0012 


0.0013 


0.0012 


0.0013 


0.0001 


0.0018 


Tagging ((Jtag) 


0.0007 


0.0016 


0.0011 


0.0014 


0.0006 


0.0012 


Background (o"bkg) 


0.0023 


0.0029 


0.0020 


0.0029 


0.0042 


0.0070 


Dependence on m^^o (ciAdep) 










0.0018 


0.0047 


Total (cJtot) 


0.0069 


0.0114 


0.0099 


0.0150 


0.0100 


0.0193 



The impact of a possible systematic mis-measurement of At (uAt) has been estimated by com- 
paring different parameterizations of the resolution function, varying the position of the beam spot, 
and the absolute z scale within their uncertainties, and loosening and tightening the quality criteria 
on the reconstructed vertex. We also estimate the impact of the uncertainties on the alignment 
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of the silicon vertex tracker (SVT) by repeating the measurement using simulated events, inten- 
tionally misaligning the SVT in the simulation. As the systematic uncertainty of the fit (cfit), 
we quote the upper limit on the bias on the and parameters, as estimated from samples of 
fully-simulated events. The model error (cTmod) contains the uncertainty on the lifetime and 
Amd, varied by the uncertainties on the world averages [HJ and also by floating them in the fit. The 
tagging error (o"tag) is estimated considering possible differences in tagging efficiency between 
and B^ and allowing for different At resolutions for correctly and incorrectly tagged events. We 
also account for uncertainties in the background (<7bkg) by varying the effective lifetimes, dilutions, 
rriES shape parameters, signal fractions, and background CP asymmetry. For the B — > Dp decay 
we also include the iiia,xiniuni bias of the a and. ci^p parameters due to the possible dependence of 
A on the vrvr'^ invariant mass (ciAdep)) as discussed in section |31 

5 PHYSICS RESULTS 

From the unbinned maximum likelihood fit we obtain the result: 

^Dtt ^ -0.013 ± 0.022 (stat.) ±0.007 (syst.) , cg^ = -0.043 ± 0.042 (stat.) ± 0.011 (syst.) , 
^D'n ^ -0.043 ±0.023 (stat.) ±0.010 (syst.) , cg*'^ = 0.047 ± 0.042 (stat.) ± 0.015 (syst.) , 
a^P = -0.024 ±0.031 (stat.) ±0.010 (syst.) , c^^ = -0.098 ± 0.055 (stat.) ± 0.019 (syst.) . 

The largest correlation with any linear combination of other fit parameters is about 20% and 30% 
for the a^" and the c^^ parameters respectively. Figures |2l to 0] show the At distribution for events 
tagged with leptons (which have the lowest mistag probability for the B^D^n^ , B^D*^tt^ , 
B^D^p^ modes). We combine our result with the result obtained on partially reconstructed 
B L»*±7r=F decays [Zj and we use the frequentistic method described in Ref. 7 to set a constraint 
on I sin(2/3 ± 7)|. The confidence level as a function of | sin(2/3 ± 7)| is shown in Fig. |S1 We set the 
lower limits | sin(2/3 ±7)] > 0.64 (0.42) @ 68% (90%) confidence level. 

6 SUMMARY 

We have studied the time evolution of fully reconstructed B'^^D^*'>^'k^ and B^^D^ p"^ decays 
in a data sample of 232 million T{AS) BB decays. CP violation arising from the interference 
of the CKM-suppressed and the CKM-favored amplitudes is expected to be small but sensitive to 
sin(2/3±7). 

The measured CP- violating parameters defined in Eq. El are shown in Section El No significant 
CP asymmetry is observed thus far. Using a frequentistic approach and combining our result with 
the BABAR result of Ref. 0, we set the limits: | sin(2/3±7)| > 0.64 (0.42) @ 68% (90%) confidence 
level. 
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Figure 2: At distribution for data having the decay mode i? D^tt^ in the lepton tagging 
category for, from upper left going clockwise, Btag = and Brec = D^tt^, Btag = B^ and 
Brec = -D+TT", Btag = B^ and Bj-ec = D~TT^, Btag = B^ and Brec = D^TT^ . The result of the fit is 
superimposed. In each plot, the dashed curve represents the background contribution. 
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Figure 5: Frequentistic confidence level as a function of | sin (2/3 + 7)|, obtained when combining 
our result with the result previously obtained on partially reconstructed B — > D*^7r^ decays [7]. 
The horizontal lines show, from top to bottom, the 68% and 90% confidence level. 
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